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Abstract 

We investigate the quantum origin of the primordial cosmological gravitational waves in 
the Brans-Dicke theory. We compute the number of gravitons Nk produced during inflation 
and the observables: power spectrum Pr, spectral index nj and energy density £lk- By 
comparison with General Relativity we see that the results for both theories are the same 
for the case of the particles number Nk and for the case of the power spectrum Pj and the 
energy density Q. only when ft) < 10. For the spectral index nj we found that when the 
Brans-Dicke coupling parameter ft) is bigger than unity the spectral index approximates the 
expression found in the General Relativity. This may awake us to the possibility that ft) 
varies with the cosmological scale. 

KEYWORDS: scalar-tensor gravity, gravitational waves 

PACS numbers: 04.30.-w, 98.80.-K 

1 Introduction 

The gravitational waves are tensorial fluctuations in the metric of spacetime. This particular 
perturbation is not explicitly coupled with the energy density and pressure of the matter of the 
Universe and does not contribute to the gravitational instability that generates the cosmologi- 
cal structures we see today. On the other side his study is of great interest because it supplies 
the specific signature of the metric theory of gravity. Nevertheless, they must have left special 
signatures in the polarization of the CMB anisotropics |T]-[3j. Moreover, since the gravitational 
waves were generated in the early Universe, before the time of last scattering of CMB photons, 
they can be a window to the primordial phase of evolution of the Universe. These waves are 
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predicted by Einstein's theory of General Relativity, but they still have to be directly detected. 
Great efforts are been done in this sense and there is hope that a new generation of experi- 
ments in space may allow this detection in some years p). The gravitational waves spectrum 
frequency extends over a wide range of interest, from 10 -18 //z to 10 8 //^| depending of the 
sources that generate those waves. Many works have been made in order to identify specific 
signatures of cosmological models in the spectra of gravitational waves, for example, in the 
case of quintessence model (5J and string cosmology |6]-[8j. 

The existence of a classical scalar field in the nature has been considered in many theories 
of gravitation that present alternatives to General Relativity. The prototype of scalar-tensor 
theories is the Brans-Dicke theory [|9- 12 1, whose Lagrangian in the Jordan frame is given by 



16n, 



(1) 



where is a scalar field that couples to gravity through the parameter ft), called Brans-Dicke 
parameter and J£ mat is the matter Lagrangian. 

The Brans-Dicke theory, as others scalar-tensor theories, can be formulated in the Einstein 
or in the Jordan frames, which are conformally related. Which conformal frame is the physical 
one is a very contentious issue p"3| . The most common argument against the Jordan frame found 
in the literature is that the scalar field energy density can assume negative values. This happens 
because of the nonmininal coupling between the scalar field and geometry. The terms with 
the second covariant derivatives of the scalar field in the Brans-Dicke field equations contain 
the connection r" v and therefore part of the dynamical description of gravity. But the energy 
density can be made nonnegative with the help of a new connection f " v [[Til which is also 
associated with the physical metric g^ v . Based on the considerations presented in the previous 
reference we work here considering Jordan frame as the physical frame. 

One can expect that the presence of the scalar field leads to different predictions about the 
results obtained in General Relativity. The standard cosmological scenario given by General 
Relativity shows some spectacular successes but also have some puzzles, such as the flatness 
problem, the horizon problem and structure formation problems p"5J . 

The value of the parameter CO can be limited by local physics experiments like Cassini 
experiment [T6}T7J: it should be huge today, the order of CO = 40000 to satisfy the experimental 
tests of General Relativity. This result reduced the interest in the Brans-Dicke theory because 
when co — > °° in the motion equations we obtain the General Relativity. The situation has 
changed with the proposal of the extended inflation [18 19 1. The idea is to take advantage 
of the time dependence of gravitational constant G to solve the bubble nucleation problem 
that arises in tradicional inflation. In the de Sitter phase, Brans-Dicke theory predicts power- 
law inflation instead of exponential; in fact, following JT8J , at the beginning of inflation the 
Brans-Dicke solutions (for p = — p) approaches the Einstein-de Sitter solution. In the second 
stage of the inflation, both the scalar field and the scale factor grow by power law rather than 
exponential. This feature prevents the so-called graceful exit problem. However, in order to 
work, the parameter must be CO ~ 24, contradicting local observations. This drawback can be 
overcome through a generalization of the original Brans-Dicke theory, allowing the parameter 
ft) to be a function of the field itself. 



lr The wave number interval is from 10 26 m 1 to 10 16 m 
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In 1999, the SN la observations |20| showed that the universe is currently undergoing accel- 
erated expansion. A possible theoretical explanation for this acceleration is the vaccum energy 
with negative pressure, called dark energy (that violates the strong energy condition). There are 
many dark energy models, for example quintessence pTfl , k-essence [|22j, phantom p3~| and 



Chaplygin gas [24|. The simplest candidate for the dark energy is the cosmological constant 
A whose equation of state is p = —p. But there are some problems with this alternative: why 
the cosmological constant is so small, nonzero, and comparable to the critical density at the 
present? So, it is interesting to study dynamical dark energy models like the Brans-Dicke the- 
ory of gravity as a possible theory of k-essence coupled to gravity since it involves the simplest 
form of non- linear kinetic term for the Brans-Dicke scalar field <$>(t) J25| . 

This revival of the Brans-Dicke theory leads us to ask if the generation of gravity waves 
can be modified through the introduction of this scalar field. Here, we propose to study the 
evolution of tensorial fluctuation in the traditional Brans-Dicke theory and its quantum origin 
and calculate the observables that can be compared with the observational data. The main 
interest in the presence of the scalar field is the possibility that it can change the evolution of 
gravity waves in comparison with the results obtained with General Relativity [26|. 

This paper is organized as follows. In section II we show the field equations, the equations 
of motion of the Brans-Dicke theory and the background solutions. In section III we obtain 
the perturbed equations of gravity waves and their solutions in terms of Hankel functions in the 
primordial phase of the universe. In section IV the quantum formulation of the gravitational 
waves in the Brans-Dicke theory is obtained. In section V we calculate the observables of the 
Brans-Dicke theory. Finally, in section VI we discuss the results and show our conclusions. 

We use in this article standard tensor notation: the greek indices run from zero to three; 
the latin indices run from one to three; the space-time metric is taken to have signature +2, or 
(— , +, +, +); covariant derivatives are indicated by semicolons; partial derivatives are indicated 
by commas and the scalar field is a time function. 

2 Background solutions 

We assume the background Universe is spatially flat, homogeneous and isotropic, i.e., it is 
described by the Robertson- Walker metric 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) , (2) 

where the a(t) is the scale factor of the Universe, and c = 1. The energy-momentum tensor of 
the background matter takes a perfect fluid form 

T^ = (p+p)u^u v +pg^ , (3) 

with an equation of state p = ap. In these expressions p is the matter density, p is the pressure 
and a is a constant. From the action ([T|) we obtain the field equations 

Rpv-^gnvR = y^v + ^^v- ^gixvM^j + ^U>^-v-g^n^j , (4) 

□0 = t-i^T . (5) 
3 + 2(0 
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Inserting the metric Q, into equations ([4]) and ([5]) we obtain the equations of motion 

„a 871; / 2 + co + 3a + 3a(o\ 6 2 & 
- 3 a = J P { 3T2^ J + V + * ' <6> 

The above equations must be supplemented by a conservation equation of the energy- 
momentum tensor 

7# v = . (8) 
When /i = we obtain the energy conservation 

p+3%(l + a)=0 , (9) 
while, for fJL = j we get the condition 

P, ; = . (10) 



Background solutions [27 1 can be obtained assuming that the scale factor a(t) and the scalar 
field <j)(t) have a power-law form 

a{t)ocf , 0(f) «J S . (11) 
By direct substitution of the above relation in the equations ([6]), ([7]), ([9]) and ( 10), we have 

-3r(r-l) = 8^pf 2 - s f \+s 2 {l + co)-s , (12) 

8tt 

,2 + , (3r _ 1) = ^_^ p? 2-, (1 _ 3a) (13) 

3 + 2(0 

p+ 3^p(l + a) = , (14) 
P,i = . (15) 

In what follows we shall consider the inflation case (a = — 1), the radiation phase (a = 1/3) 
and the dust era (a = 0). In the first scenario we have a drastic expansion of the Universe during 
the early period of the Big Bang. In this case, the background solutions are 

s = 2 , r =co + ^ . (16) 

In the radiation phase, with a = 1/3, we have the same background solution that the General 
Relativity 

* = , r= l - . (17) 

In the dust era (a = 0) we obtain 

1 + 0) 

s = 2-3r , r = 2-—— , (18) 
4 + 3(0 
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3 The gravitational waves equation in the Brans-Dicke the- 
ory 

Cosmological gravitational waves are obtained by means of a small correction hij in equation 
(|2]), which represents the metric. Hence, the general expression of the metric ([2]), related to the 
unperturbed metric, is replaced by 



cb 2 = -dt 2 + [a 2 (t)8ij + hjj] dx'dx j . 
Appling the metric ( 19) in the field equation ([4]) the resulting expression is 



<j) d\ . 
- - - hij + 
<p a J 



k 2 „a „d6 

2 2— 

a 2 a aty 



hij = 



(19) 



(20) 



where k is the wavenumber, the dots indicate cosmic time derivatives and we have written 
hjj(t,x) = h(t)Qjj, where Qy are the eigenmodes of the Laplacian operator, such that Qa = 
Qki,k = 0. 

Performing the transformation from the cosmic time to the conformal time, a(t) dr\ = dt, 
and representing the derivatives with respect to r\ by primes, equation ([20]) assumes the form 



K'j+ 



a 



(j) a 



h'ij + 



, 2 a" a' 2 a'(j)' 
k 2 -2— + 2^-2—^ 
a(p 



hij = 



a or 

To resolve this differential equation we need to perform the Fourier transformation 



hij(x,ri) = V\6n £ 

A=(g>,6 



d3k ^) (fr A^H^) r -ilt 



(21) 



(22) 



where the polarization tensor (k) can be decomposed as (k)s lJ '^(p) = 28^' 8^ (k — 
p) to the two polarization states <8> e ©. Then, we have 



■ k2+ a^ + l_r_l(l\ 2 + aW~ ] 
a 2 (j) 4 J a(j) 



M(A)(n) = 



(23) 



In this way, we rewrote the equation (21 ) in terms of a simple harmonic oscillator. 

We can express the equation (23) in terms of the new parameter 4>(t]) = a(r\)yjty(r\). 
Hence, the equation ([23]) becomes 



M(A)(n) = o • 



With the background solutions of the scale factor and the scalar field we get 



(24) 



<E"(t]) 



12 » 



G(co) 



if a = — 1 with — °° < 71 < — 771 
if a = 1/3 with -T]i < 77 < 772 



(rt,F . if a = with t] > Th 



where 



and 



r ldust — 12 



F(co) 
G((o) = 

2 + 2(0 



2(3 + 2©) (1 + 2©) 

(l-2(0) 2 

(3 + 2(0)(l + (0) 

(2 + (0) 2 ' 



1 ©-1/2; V m 



2 + (0 

The solutions, for each stage of evolution of the Universe are 
1. Inflation era (a = — 1) 



(25) 



(26) 



^(kr 1 ) = ^n[A v jr v {1 >(k\r,\)+ A 2 J?r> (k\r,\) 



(2), 



(27) 



where A\ and A2 are integration constants, J^f^ and ^f( 2 ) are the Hankel functions of 
first and second kind, respectively, and v = ^F(co) + 1/4 is the order of the Hankel 
functions. 



2. Radiation era (a = 1/3) 

ll(kri) =B ie ik7] +B 2 e- ik ^ , 
where B\ and B 2 are yet integration constants. 

3. Dust era (a = 0) 

H(kr)) = y/rj C h ^ l) (k(r, - f\ dust )) + C 2 m£ 2) [k{r] - f\ dust )) 



(28) 



(29) 



where C\ and C 2 are integration constants and a = a/G((o) + 1/4 is the order of the 
Hankel functions. 



4 The quantum gravitational waves equation in the Brans- 
Dicke theory 



In this process we start with the Lagrangian density of the gravitational waves in the Brans- 
Dicke theory, given by 

eh ..... . .. 

(30) 



that leads to two corresponding Hamiltonians 



» m - -If* 

" m - -y 



X 



d 3 * 



k z -2—h7C + 8 ij HjLIj 



(3D 
(32) 
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where % and % are the canonical momentum. The above Hamiltonians differ each other by 
addiction or subtraction of a total conformal time derivative. 



The expression of the quantum Hamiltonian obtained by the equation (32) for the gravi- 
tational waves in the Brans-Dicke theory is constructed by using the canonical commutation 
relations 

[fl(x 1 r 1 ) 1 7t(y 1 r 1 )] = i8^(x-y) . (33) 

So, by transforming the normal modes of oscillation of the gravitational waves in field 
and momentum operators in the Heisenberg description of quantum mechanics we obtain the 
Hamiltonian 

tf(2)^#(2) = _. 



d 3 jc 



(34) 



where the field operator is constructed from its classical field as follows 



htj->fiij(x,ri) 



Vl6n 



1 



£^\k)d 3 k 

(2^)3/2 



k,{X) 



(35) 



with the operator given by 

A| = %(t?o)A(t?)+«^(t]o)/,*(t?) 
and the function satisfying the following equation 



fk+ k 



fk = 



(36) 



(37) 



that is the same differential equation that describe the behavior of the classical function jU(tj). 

Now, to obtain coherent states of the relic gravitons in the Brans-Dicke theory we consider 
the creation and annihilation operators 



'k I A i 



(38) 



whose temporal evolution is given by the Heisenberg equation and by the Hamiltonian (31 ), 
writen now as operator 



(i) 



d 3 x 



(39) 



with the general solutions given by 

= u k {T))a l {7] Q )+ v fc (77)a t 7 (Tjo) 



fl^(T]) 



V*(T?)^(no) + "i(T?)^(TJo) 



where T]o is a fixed initial time and the functions m^(t]) and v k (r\) behave as 



du k <t>\ 
-— = iku k + —v k 
dr\ 4> fc 



— =ikv k + —u k 
drj 4> 



(40) 



(41) 



The equations ( J40| ) are the Bogoliubov transformations [28j of the gravitational waves in 
the Brans-Dicke theory and u^rj) and v^rj) are the Bogoliubov coefficients that satisfy the 
relation 



\Uk\ 



\ v k\ 



1 



(42) 



result that guarantees the unity of the temporal evolution of these operators. 

The number operator is defined as N k = v\ v k- So, in our model the number of gravitons is 
given by 



N t 



+ 



i 

2k 



dfk 



dr} 



1 3>' 
2k® 



dr] 



k +fk 



dfk 
drj 



1 

2 



(43) 



where the relation between v k an d fk is 



V2kdri 

When krj > > 1 the number operator simplifies as 



k — i 
2k V 4> 



(44) 



N k nk\f k \ 2 1 + 



1 

2k 2 V 4> 




1 

2 



If ^7] << 1 we have 



N k *k\f k \ 



1 

2 



(45) 



(46) 



In Fig. [T]- Fig. |4]we show the behavior of the number of particles in the inflation era in the 
Brans-Dicke theory to different values of CO and small and big values of the wave number and 
we compare with the General Relativity. The aim here is, first of all, see how the Brans-Dicke 
model behaves for different values of 0), since local observations show that CO is big and, as we 
have written before, small values should be more interesting in large scales. Secondly, when 
gravitational waves of cosmological origin are observed, the Brans-Dicke theory and General 
Relativity may be confronted regarding the evolution of the tensor modes. 



5 Observables in the Brans-Dicke theory 

In this section we apply the formulation of quantum gravitational waves in Brans-Dicke theory 
to determine the most common and interesting observables in cosmology. The experimental 
data at present do not allow us to include or exclude the presence of a primordial spectrum of 
relic gravitons compatible with any model. The idea is that these theoretical results calculated 
here might in the near future be compared with the observational data obtained by the projects 
related to the search for gravitational waves of cosmological origin. 
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Figure 1: The number of gravitons produced by inflation in the Brans-Dicke theory for different 
values of ft), and also in the General Relativity, for the wavenumbers k = 10~ 5 rrC 1 and k = 
10" 1 mT x . 
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Figure 2: Same as the previous figure, but for wavenumbers k = 10 6 m 1 e k = 10 7 m 1 . 
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Figure 3: The number of gravitons for small values of ft) at two different comoving times. 



5.1 Power Spectrum 

The power spectrum is related by the Fourier transform of the two-point correlation function of 



the tensor modes in the following way [29 1 



(47) 




Figure 4: The number of gravitons for large values of CO and for General Relativity. 



where Pr{k, T]) is the power spectrum and the state |0) is annihilated by a-^ . 

The process of quantization of gravitational waves makes the classical field /? ;; - to be trans- 
formed into an operator This operator is given by the equation (35 ) where the vacuum state 
|0), in which there are no particles, is annihilated by the annihilation operator 

<%|0)=0 . (48) 

The calculations of the expectation value of the field operators yields 

d 3 * lA(n)l 2 



(o| hi^nw'wm = 64»y w ^ 



-ik- {x—y) 



(49) 



Knowing that gravitons can be produced in isotropic models and using spherical coordi- 
nates, we have 



(0\hl(x,ri)h i j(y,ri)\0)=647i 



k 2 \f k \ 2 2sm(kr) 
(lit) 2 a 2 § kr 



dk , 



(50) 



Comparing the equation (47) with the previous equation we have, within the Brans-Dicke 
model, the power spectrum of the tensor modes parametrized as 



k 3 I ft I 2 

JV(M)=32jA^ 
In the inflation phase the function f k (k\r\\) is defined as 



(51) 



fk{k\n\) 



K (k\7}\) 



(52) 



and by using the small argument limit of the Hankel functions |30|, we have 

8tz;2 2v T 2 (v) k 2 



mm; 



a 2 (j) 



(k\v\) 



l-2v 



(53) 
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Figure 5: The power spectrum of the gravitational waves in the Brans-Dicke theory as a function 
of wavenumber for small values of ft). 
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Figure 6: The power spectrum of the gravitational waves in the Brans-Dicke theory for large 
values of ft). 



By definition H = - = X^, and the conformal time in this model can be written as 



I 7 ? I 



2(0 + 1\ 1 



2ft) - 1 / aH 



(54) 



Thus, the power spectrum in the Brans-Dicke theory during inflation for the modes out of the 
horizon becomes 



mm; 



8tt2 2v T 2 (v) {2(D + l\ l - 2v H 2 ( k x 3 - 2v 



71' 



2ft) -1 



(j) \aH 



(55) 



In Fig. [5]- Fig. [6] we show the behavior of the power spectrum Pr{k, r\) as a function of the 
wavenumber and of the parameter ft). 



5.2 Spectral Index 

The spectral index is obtained by the expression 

dlnP r (M; 



nj 



dlnfc 
11 



(56) 



where a = aH = ( §§jrj- ) ^ * s tne moment when the tensorial modes cross the horizon again. For 
this calculation we introduce the parameter £, that measures the rate of decrease of the Hubble 
parameter during the inflation phase and it is given by 



£ 



1 



dt \H 

In the Brans-Dicke theory we obtain that 



1 dH 
a~H^dr\ 



(57) 



2(0 + 1 



(58) 



Thus, we have 



dln£ 



In 







t kdH dln<A 
'H~d~k~ dlnfc J 



(59) 



The first term on the right side is calculated as follows 



kdH 
'H~dk 



k dH d\r]\ 



Hd\r}\ dk 



k dH fd\r] 



and, by using the equation (54), we have 1 3 1 1 



d|rj| 



dk 



d 
dk 



2co+l\ 1 



2(0-1 J aH 



H drj \ drj 



1 

dk \k 



d\rj\ 



dk 



(60) 



1 

1? 



(61) 



In this way, with the help of the equations ( |58| ), ( [60] ) and ( |6T[ ), the equation p9\ can be 
written as 



n T = -2e 



2(0+1 
2(0-1 



dln0 



dink 



(62) 



When (0 — > °° and is constant we obtain the same result that in General Relativity [31 1, i. 
e., nj = — 2e. 

Likewise, we can calculate the term dependent on the scalar field in equation ( [62] ) 

-l 



dln0 



dln£ 



kd§ 




k d0 d?7 




k d0 




dk 


a 


^ dr\ dk 


a 


" <pdt] 





ci 1 77 1 



(63) 



that with the help of equation (54) and with the background solutions of the scalar field we 
obtain the following result for the spectral index nj 



-2e 



2(0 + 1 
2(0-1 



2(0-1 



(64) 



Again, the above result gives the same value of the spectral index of the tensorial modes of 
General Relativity at the time of inflation, when CO — > °°. 
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We can write the expression of nj in terms of the parameter CO. For that we replace the 



relation (58 ) in the result found above, such that 



By using (57 ) we obtain 

2 

2(0 = --1 , (66) 
e 

which can be replaced in the final result for the spectral index such that 

Ae 

n T = = -4e(l-e)" 1 . (67) 

£ — 1 

For small values of £ we can expand the term (1 — £) _1 obtaining 

n r = -4£(l+£ + ^(£ 2 )) , (68) 
and considering only zero-order terms we have 

n T -4e . (69) 

5.3 Energy Density 

With the Lagrangian density of the gravitational waves in the Brans-Dicke theory, equation ( [30] ), 
we obtain the energy-momentum tensor of the gravitatinal waves 

T ^ = "L\ (^Mv^ 7 - Ign^dahijdph^ . (70) 



Considering 



we find that 



where 



hij = Vl6n £ ,h^ = ^ , (71) 

7> = -2$ (dphdyh - yn V g aP d a hd p h\ , (72) 



/; // // . (73) 



The energy density is the (00) component of the energy-momentum tensor of the gravita- 
tional waves ^ 

p = T ° = <l>(^+g i jd i hd j hj . (74) 
To calculate the quantum energy density we first substitute the following 

h^h = ^— , (75) 
13 



such that the expectation value of the energy density is given by 



(p) = (0| p |0) = 4 «0| h'h'* |0> + (0| Vh ■ Vh* |0» 



(76) 



With the help of the expansions 

1 r d 3 k 



A 2/(2^)V2 

^2-y 



(2^) 3 /2 



(77) 



the solutions 



and the commutation relations 



% = ^(»to)/*(»?)+^j(»?o)Jr(n) , 

% = %( 7 7o)^(r?)+a t r (r7oM(r?) , 



%4 



5< 3 >(*-*) , 



(78) 
(79) 



we obtain 

(0\h'h'*\0) 



1 /■ d 3 fc 



/\ 2 



(0|V/z-V/**|0) = 



1 /■ d 3 fc 



4/ 



a 2 <M (2^:) 3 



(81) 



Thus, the expectation value of the energy density of the gravitational waves in the Brans-Dicke 
theory can be write as 



(P) 



1 r d 3 k 



a 4 J (2k) 3 



\gk\ 2 +[k 2 + 



l\ 2 N 



2 & 

\fk\ --^{gkfk+Skfk) 



-ik(x—y) 



(82) 



Considering the presence of an isotropic background of relic gravitons we obtain 



If k 2 dk 
{p) = ^Jj2nJ 3 



l**r+U 2 + 



/\ 2 N 



2 4/ * 
- ^ (Skfk + 8k fk) 



2sinfcr 
kr 



(83) 



where r = \x — y\ . In the limit x — > y we have 



1 /-fc 2 dfc 



2^2 



l**f+ r + 



/\ 2 N 



(84) 



The energy density of the relic gravitons per logarithmic interval of the wavenumber will be 
given by 



d(p) k 3 



dink 2n 2 a A 



\gk\ z +\k z + 



\fk\ Z -^(gkfk*+fkg*k) 



(85) 
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The spectral energy density per logarithmic interval of the wavenumber is defined as 

p c din A: 



i. e., 



g(m; 



3flJ 2^ 2 a 4 



»r + l^+(*J \fk\ --^(gkfk+fkgk) 



(87) 



where p c is the critical density of the Universe. 

For all the oscillation modes that are inside the Hubble radius (i. e. krj » 1), we have 
gk « Tikfk- This implies that 

and, rewriting in terms of the power spectrum we have 

k 2 _ / 1 /4>^ 2N 



Q(M)w i^M 1 + 2*UJ J ■ (89) 

ct>' 

For modes outside the Hubble radius (i. e. kr\ « 1) we have gk~— i^fk- This result allow 



us to find that 



Snk 5 , ,. 2 



which can also be rewritten as a function of the power spectrum P T 

aM ~2^p PT(k ^- <91) 

We stress that when the modes re-enter the horizon the energy density is proportional to the 
power spectrum. The behavior of the spectral energy density are shown in Fig. [7]- Fig.[9j 



6 Conclusions 

In this work we have considered the quantum gravitational waves production during inflation in 
the Brans-Dicke theory. In general, we realize that, in contrast with the General Relativity, not 
only the scale factor a(t) but also the Brans-Dicke scalar field <j>(t) models the time evolution 
of the gravitational waves of cosmological origin. 

In the initial inflationary phase of the Universe there is a natural initial vacuum state, from 
wich the particle occupation number can be determined. We computed the quantum particles 
number created during extended inflation and the results indicate that the number predicted 
by this model approximates that predicted by General Relativity when the Brans-Dicke param- 
eter CO is big enough, as we can see in Fig.[TJ- Fig.|4j 
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The graviton energy density in the Brans-Dicke theory as a function of k, for small 
co and fixed r\ . 



1,1=10- 




5x10' 10 5 
k (m _1 ) 



1.5x10 s 2xl0 5 



[J7|=10- 




5x10- 



10 + 1.5xl0 4 
k (in -1 ) 



2xl0 4 



Figure 8: The energy density of the gravitational waves in the Brans-Dicke theory as a function 
of k compared with the energy density in the General Relativity, for some different values of CO 
and fixed r\ . 
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Figure 9: The energy density of the gravitational waves in the Brans-Dicke theory as a function 
of | T) | compared with the energy density in the General Relativity, for some different values of 
CO and with fixed wavenunber k. 



Afterwards we calculate three observable parameters: the quantum power spectrum P^, the 
quantum spectral index nj and the quantum energy density £l{k,r\) of the gravitational waves. 
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The result obtained for the power spectrum shows that, different from the local observations 
and from the previous result acquired for the graviton number N^, this model is comparable to 
the General Relativity one only if co assumes small values like 0.5 < CO < 1. Big values of CO 
make the power spectrum to explode. 

For the spectral index we obtain a different result. We can compare the outcome of the 
General Relativity for the spectral index nj with that of the Brans-Dicke theory when CO » 1, 
result that agree with the one obtained for the quantum particle number N^, and when co — > °° it 
yields a scale invariant perturbation. 

In the case of the quatum energy density big values of CO make Q.(k,r}) to explode and the 
General Relativity prediction is approximated when 1 < CO < 10, depending on the values of 
the wavenumber and the conformal time. 

We think the main problem with the scenario of Brans-Dicke is valued at CO. In our study, 
the values vary significantly from the values obtained by traditional testing sites. The param- 
eter co assumes small or large values depending on the observable calculated. However, these 
drawbacks can be, in principle, circumvented if the gravitational coupling is scale-dependent, 



as suggested by considering quantum effects [32 33 1. In this work, we made no attempt to rec 



oncile the values of CO deducted from local tests, with the corresponding values deduced from 
the analysis performed here. But, evidently, the problem must be addressed in order to have 
a completely realistic cosmological scenario. When gravitational waves are detected we can 
verify if CO has different values for different cosmological scales. 

A complete study of the Brans-Dicke theory in terms of fluctuations of quantum origin also 
involves the calculation of scalar perturbations. This is a natural continuation of this work. The 
results of scalar perturbations during inflation using the Brans-Dicke theory can be compared 
with the results of General Relativity and with the available observational data, especially data 
from the cosmic microwave background radiation. We can reduce the number of the variables 
calculating the ratio between the tensor and the scalar power spectra rj = Pj/Ps and with this 
results to obtain direct information about the state of the inflationary universe. An other step 
to generalize our study is to do the same calculating for both scalar and tensor perturbation but 
in a model where the Brans-Dicke parameter CO varies with the scalar field 0, i. e., CO = Co((j)). 
Such work is in preparation and will be available soon. 
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